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In this paper we extend and unify the results of |£C| and | |l9| . As a consequence, the results of 
|2o[ are generalized from the framework of ideal polyhedra in H"^ to that of singular Euclidean 
structures on surfaces, possibly with an infinite number of singularities (by contrast, the results 
of pol can be viewed as applying to the case of non-singular structures on the disk, with a finite 
number of distinguished points). This leads to a fairly complete unHerstandirig nf t he m oHuli space 
of such Euclidean structures and thus, by the results of Epstein, Penner, {||lq; pt h_3|) the author 
|l9| Q , and others, further insights into the geometry and topology of the Riemann moduli space. 

The basic objects studied are the canonical Delaunay triangulations associated to the afore- 
mentioned Euclidean structures. 

The basic tools, in addition to the results of |19| and combinatorial geometry are methods of 
combinatorial optimization - linear programming and network flow analysis; hence the results 
mentioned above are not only effective but also efficient. Some applications of these methods to 
three-dimensional topology are also given (to prove a result of Casson's). 
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Introduction 

In the paper we gave the following description of the angles of ideal polyhedra 
in H'^: let P be a combinatorial polyhedron, and let A : E{P) [0,7r) be a 
function. Then, there exists an ideal polyhedron combinatorially equivalent to P, 
such that the exterior angle at every edge e is given by A{e) if and only if the sum 
of A{e) over all edges adjacent to a vertex of P is equal to 2tt, while the sum of 
A{e) over any nontrivial cutset of edges of P (that is, a collection of edges which 
separates the 1-skeleton of P, but which are not all adjacent to the same vertex) is 
strictly greater than 27r. Furthermore, it was shown in [|l9| that the dihedral angles 
determine the ideal polyhedron up to congruence. 

It was observed in |18; |l^ that this was a special case of the problem of char- 



acterizing the Delaunay tessellations of singular Euclidean surfaces - there is a 
canonical way to associate ideal polyhedra to Delaunay triangulations of a convex 
flat disk with convex polygonal boundary. The general situation is described in 
detail below, but one of the goals of this paper is to extend the characterization 



^This is a slightly modified version of a preprint first distributed in September 1996 
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above to the completely general case of singular Euclidean surfaces with boundary. 



This is duly done, see, eg. Theorems 3.1, 4.8, 4.11 



Consider a surface S, equipped with a Euclidean (or, more generally, a similarity) 
structure E, possibly with cone singularities. Assume that there is a discrete collec- 
tion P = {pi, . . . of distinguished points on S, and assume that P contains the 
cone points of S. There is a canonical tessellation attached to the triple {S, E, P) 
- the so-called Delaunay tessellation (see, eg, |l^). The moduh space M. of such 
triples is then naturally decomposed into disjoint subsets M.t, corresponding to the 
different combinatorial types T of the Delaunay tessellation. This is a canonical 
decomposition of A^. In the paper ||l^ I studied the subsets A^Ti and showed that 
the dihedral angles of the Delaunay triangulation are natural coordinates (moduli) 
for A^T, which induce on AAt the structure of a convex polytope. 

The aformentioned decompositon of moduli space then becomes a polyhedral 
complex, the top-dimensional cells of which corresponds to Delaunay tessellations 
which are triangulations, while pairs of adjacent top-dimensional cells differ combi- 
natorially by a diagonal flip. This decomposition is closely related to the well-known 
Harer complex (see, eg, |l^]). As mentioned above, the top-dimensional cells of 
this complex are identified along some of their lower-dimensional faces, while other 
lower dimensional faces correspond to degenerations of the Euclidean structures 
of {S,E,P). It is then clear that the polyhedral structure of the cells A4t is of 
considerable interest. However, in [ p^ only an indirect description was given - A4t 
was shown to be a convex polytope by virtue of being an image of another convex 
polytope under a fairly complicated linear map. The methods of come from 
hyperbolic geometry and are based on the study of dihedral angles of compact hy- 
perbolic polyhedra in [ p^ , so do not easily generalize to the case of general singular 
Euclidean and similarity structures alluded to above. In the current paper, methods 
of mathematical programming and the results of |19| are used to give a completely 
general extension of the result of (described in the beginning of this Introduc- 
tion) to Delaunay triangulations of arbitrary singular surfaces (Theorems 3.1, 4.8, 
4.11). Since the arguments do not depend on the results of |2^, we have a different, 
essentially com binatorial, proof of a principal result (Theorem 0.1) of that paper 
(Theorem 4.11 here) . The other result of Ig^l ~ the characterization of finite- volume 
polyhedra - is, seemingly, not accessible by current methods. The above-mentioned 
result permit us to get good understanding of the boundary structure of the cells 
Mt, and, consequently, of M itself. Since A4t fibers (in multiple ways) over the 
moduli space of finite area hyperbolic structures on 5, some information is obtained 
about the latter moduli space. 

In addition, the methods, combined with a geometric estimate, allow us to give 
a description of dihedral angles of Delaunay tessellatio ns o f (5, E, P), where (S, P) 
is not necessarily of finite topological type (Theorem 6.1). This stops well short 
of solving the moduli problem, unlike in the finite case, but a conjectural picture 
seems fairly clear. 

The methods are also brought to bear onto some questions in combinatorial 
geometry, and to provide efficient algorithms for solving the "inverse problem" of 
determining when a combinatorial complex, or a combinatorial complex equipped 
with dihedral angle data, can be realized as the Delaunay tessellation of a singular 
Euclidean surface. 
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In addition, we use our methods to prove some observations of Casson on ideally 
triangulated 3-manifolds. That subject is not so far removed from the geometry of 
similarity and Euclidean structures on surfaces. Indeed, the basic idea of [|l9| is to 
study the similarity and Euclidean structures by means of constructing a canonical 
hyperbolic polyhedral complex as a "cone" over the surface being studied. 

The plan of the paper is as follows: In section ^ the relevant definitions and 
results of |l^] are recalled. In section |^ we describe a set of constraints which must 
be satisfied by the dihedral angles of any (not necessarily Delaunay) triangulation. 
In section ^ we show that these constraints are actually sufficient under the as- 
sumption that the triangulation is Delaunay, and refine them to a minimal set of 
constraints. In section we show how the results apply to ideal polyhedra, and 
in particular to characterize infinite ideal polyhedra. In section ^ we comment on 
the boundary structure of the moduli space of singular Euclidean structures, and 
describe a correspondence between the Euclidean and hyperbolic structures, which 
hopefully clarifies the picture. In section |l^ we apply the methods of section || to 
the study of ideal triangulations of 3-manifolds. In section |^ we give a network 
flow interpretation of the results of section ^. In addition to the intrinsic interest, 
this allows us to give efficient algorithms for deciding whether a weighed graph is 
the 1-skeleton of a Delaunay triangulation (with weights being the dihedral angles). 
These computational issues are discussed in section ||. In section ^ we discuss some 
combinatorial-geometric applications of the results of section |^. 

1. BACKGROUND 

1.1 Singular similarity structures. 

Consider an oriented surface S, possibly with boundary, and with a number of dis- 
tinguished points {pi, . . . ,p„}. A similarity structure on S is given by an atlas for 
S\{pi, . . . ,Pn}, such that the transition maps are Euclidean similarities. A similar- 
ity structure induces a holonomy representation ff^ of F = 7ri(S'\{pi, . . . ,Pn}) into 
the similarity Sim (E^) ~ C*, where the tilde indicates the universal cover. We 
define the dilatational holonomy as the induced representation Hd : F ^ R, where 
Hdil) — log dilatationi/s (7) = \og\Hs{'^)\. The rotational holonomy can almost be 
defined as Hj. (7) = arg Hg (7) , but for the slight complication that we need to take 
the argument in C*, since we want to distinguish, eg, the angle of 47r from one of 
27r. This notion of argument is what is used in the sequel. In particular, if 7 is a 
loop surrounding one of the distinguished points pi, then i?r(7) is the cone angle at 
Pi. In the case where H^iT) = {0}, the similarity structure is a singular Euclidean 
structure, with cone angles defined by Hr as above. In the sequel, references to 
the holonomy of a similarity structure will actually mean the dilatational holon- 
omy H4. A more concrete way to think of both similarity and Euclidean structure 
is by assembling our surface out of Euclidean triangles, in the pattern given by 
some complex T. In the case where the lengths of the edges of the triangles being 
glued together agree, then we have a singular Euclidean structure. If not, then we 
have a similarity structure. In either case, the vertices of T are potentially cone 
points, with cone angles given by the sums of the appropriate angles of the incident 
triangles. 
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Consider an oriented surface S, possibly with boundary, and assume that S has 
a EucHdean metric (or, more generally, a similarity structure - with the exception 
of the results of Section ^ , some of which depend on the metric Theorem 6.3, 



the metric structure or lack thereof plays no role in the arguments), with cone 
singularities. We will be dealing with semi-simplicial triangulations of S, that is, 
triangulations where two and where distinct closed cells might intersect in a collec- 
tion of lower-dimensional cells. All triangulations will be assumed semi-simplicial, 
unless specified otherwise. A subcomplex of T will be called closed if whenever 
an open face F is in T, so are all of the faces of dF. 

Assume now that the surface S is equipped with a finite geodesic semi-simplicial 
triangulation T, such that the 0-skeleton of T, which is denoted by V{t), contains 
all of the cone points of S. Each face of T is then a Euclidean triangle. There are 
two kinds of edges of T - the interior edges, incident to two faces of T, and the 
boundary edges, incident to only one face of T. 

Definition 1.1. Let e be an edge ofT. First, suppose that e is a boundary edge, 
and let t — ABC be a face ofT incident to e, so that e = AB. Then the dihedral 
angle i5(e) at e is the angle oft at the vertex C. Now, assume that e is an internal 
edge ofT, so that e is incident to ti = ABC and t2 — ABD, so that e ~ AB. Then 
the dihedral angle at e is the sum of the angle of ti at C, and the angle of t2 at D. 
The exterior dihedral angle S'(e) at e is defined to be S'{e) = tt — (5(e). 

The cone angle at an interior vertex v of T is the sum of all the angles of the 
faces of T incident to v at v; the boundary angle at a boundary angle is defined in 
the same way. 

The following is a slight extension of |jl^. Lemma 4.2]. 

Lemma 1.2. The cone angle at an interior vertex v G V{T) is equal to the sum 
of the exterior dihedral angles at the edges ofT incident to v. At a boundary vertex, 
the boundary angle is equal to the sum of the exterior dihedral angles as above, less 

TT. 

Proof. First, let v be an interior vertex. Suppose that there are n triangles 
ti, . . . ,t„ incident to v. The sum of all of their angles is rnr. The cone angle at v 
is the sum of the angles of the triangles ti at v. On the other hand, the sum of the 
dihedral angles of the edges ei , . . . , e„ incident to v is the sum of all of the angles 
of ti not incident to v. Thus, 

n 

niT = Cone angle at v + ^ S{ei). (1) 
1=1 

The result follows by rearranging the terms. 

If w is a boundary vertex, and there are n faces incident to v, then there are n+1 
edges, and equation |^ becomes: 

n+1 

rnr = boundary angle at v + S{ei), (2) 



and the result follows by rearranging terms, as above. □ 
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Observation 1.3. The sum of all of the dihedral angles of all of the edges of 
T is equal to the 7r|y(T)| - in combination with the Lemma above this gives the 
Gauss-Bonnet theorem in this polyhedral context, since the curvature at an interior 
vertex V ofT is defined to be 2-K—Cone angle at v, while the curvature at a boundary 
vertex is defined to be t: — boundary angle at v. 

Proof. Every angle of every face of T is opposite to exactly one edge of T . □ 

The next theorem is Theorem 6.16]. 

Theorem 1.4. Let A : E{T) (0, 27r) be an assignment of dihedral angles 
to the edges of T , and Hd a holonomy representation. There exists at most one 
singular similarity structure on S with holonomy Hd (and in particular, at most 
one singular Euclidean structure, up to scaling), such that so that 5{e) — A(e), for 
every edge e G E(T). 

Definition 1.5. A triangulation T with (5(e) < tt for every interior e £ E(T) 
is called a Delaunay triangulation. 

Let A be a map - A : E{T) — ^ (0, 2tt). When does there exist a singular Euclidean 
metric on S with the dihedral angles prescribed by A? It is clear that there are 
certain linear constraints which must be satisfied - to wit, for every face t = ABC 
of T, we must be able to find angles a, P, and 7, such that: 

Positivity. All angles are strictly positive. 
Euclidean Faces. For every face t, at + Pt + = tt. 

Boundary edge dihedral angles. For every boundary edge e — AB, incident to 
the triangle ABC, 7 = A(e). 

Interior edge conditions. For every interior edge e — AB, incident to triangles 
ABC and ABD, a + S = A(e). 

The system of equations and inequalities above specify a linear program L. The 
feasible region (set of solutions) of L must be non-empty in order for us to have any 
hope of having a singular Euclidean metric on S with prescribed dihedral angles. 



One of the principal results (Theorem 6.1) of |19| is that if all of the dihedral angles 
are no greater than tt (that is, the triangulation is Delaunay), then Conditions 1-4 
are also sufficient, thus: 

Theorem 1.6. If the feasible region of L is non-empty and every dihedral angle 
is at most tt, then there exists a similarity structure with any prescribed holonomy 
Hd and, in particular, a singular Euclidean metric on S with the prescribed dihedral 



angles - this structure is unique by Theorem L4 (and the metric is unique up to 
scaling). 

2. NECESSARY CONDITIONS ON DIHEDRAL ANGLES 

In order for the linear program L to have any chance of having a solution, the 
dihedral angles A must satisfy some constraints. Indeed, suppose L has a solution. 



Condition 2.1. All of the dihedral angles must be positive. 
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Furthermore, the sum of aU of the dihedral angles of all of the edges of T must be 
equal to the sum of all of the angles of all of the faces of T, or F{T)-k. On the other 
hand, it is almost equally obvious that if ii, . . . , t„ e F{T) is is some proper subset 
of the faces of T, then the sum of the dihedral angles at the edges incident to one 
of the U must be strictly greater than the mi. In other words: 

Condition 2.2. Let T C F(T), and let E{T) he the set of all edges incident to 
an element of T . Then 

E ^(e)>'^l-^l' (3) 

with equality if and only if T — F{T) or E{J-) ~ 0. 

Definition 2.3. For a subcomplex T ofT, we define the excess of J-' to be 
excess !F= E ^(e)— 7r|jr|, 

e6_E(JP) 

Definition 2.4. An edge e of a subcomplex J- of T is a relative boundary edge 
of T if it is incident to a top- dimensional face of !F on exactly one side, and is not 
a boundary edge of T . 

The following lemma will prove useful in the sequel. 

Lemma 2.5. Let A : E(T) — > (0,7r] be an assignment of dihedral angles to edges 



of T satisfying condition 2.L Let ti,t2, ■ ■ ■ ,tn be a collection of closed triangles of 
T , and let = ti U t2 D . . . U tn - assume T ^ T . Let 'Y^q T be the sum of the 
dihedral angles of the boundary edges of T . Then 

< excess T <^T . (4) 



Proof. The first inequality of ^ is a restatement of Condition 2.2, applied to 
the subcomplex T . To show the second inequality, let 

U ^ 

tGF{T)\{tt,....t„} 

Evidently, 

J^uT = T, 

while 

Tr\T=dT. 



Applying the conditions 2.2 to JT, we see that 



0< E A(e)-7r|F(^)|, (5) 

while applying them to JT, we see that 

0< E A(e)-7r|F(:F)|. (6) 



7 



Note now that 

\FiT)\ + \F{T)\^\FiT)\, (7) 

while 

A(e)+ J2 Me)-Y.^= E ^(e)- (8) 

Adding inequahties |5| and H, and applying equations and |[ we obtain 
A(e)-7r|F(^)| 

eef;(:p) 

< I A(e)-^|F(^)|) + I Y He)-AF{y)\ 

Y A(e)+ Y ^{e)\-{AF{T)\+^\F{y)\) 
eeE(r) eeE{T) / 

= Y A(e)+5]^ + 7r|F(r)|=5]^, 

eeE(T) d d 

where the last equality is obtained by applying the equality case of Condition 
□ 

3. SUFFICIENCY OF CONDITIONS 2.1 AND 2.2 



2.2 



Somewhat surprisingly, the trivial conditions 2.1 and 2.2 guarantees that the linear 
program L has a solution: 

Theorem 3.1. Let T be a semi-simplicial triangulation of a surface S, possibly 
with boundary, let A : E{T) (0, tt] he a function on the edges of T, a rep- 
resentation '!Ti{S\V{T)) — )■ R. There exists a similarity structure with holonomy 
Hd on S , with the Delaunay triangulation of {S, £) combinatorially equivalent to 



T , and with dihedral angles given by A if and only if conditions \2.:\ and 2.i are 
satisfied. 



We will use the Duality Theorem of linear programming to prove Theorem 3.1 
Before stating the Duality Theorem we need to recall some notions: A linear pro- 
gram L consists of a collection C{L) of constraints^ which are linear equations or 
(nonstrict) inequalities, and the objective function F{L), which is a linear function. 
The set of points in R" satisfying the constraints C(L) is called the feasible region 
of i, which is, by definition, a polyhedral region, possibly unbounded, possibly not 
of full dimension, and possibly empty. The solution of the linear program is a point 
where the objective function attains an extremum (which may be a maximum or a 
minimum). The value of the objective function at the solution is the objective of 
L. If the feasible region of L is empty, the program is said to be infeasible. Now 
we can state the Duality Theorem: 

Theorem 3.2 (Duality Theorem of Linear Programming). Let F be a 
linear program of the form: 
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Minimize c^x, subject to the constraints Ax — a, x > 0. 

Then the dual of P is the program P* : 

Maximize a-^X, subject to the constraints A"*" A < c. 

The feasible region of P is nonempty if and only if the objective of P* is bounded. 
Conversely, the feasible region of P* is nonempty if and only if the objective of P 
is bounded. If neither feasible region is empty, then the values of the objective 
functions of P and P* are equal. 



Remark. While the primal program P in the statement of Theorem |3.2| appears 
to be of a very special form, it is not hard to see that any linear program can 
be written in this form. Indeed, if our linear program asked us to maximize the 
objective, we can always convert it to a minimization problem by multiplying c 
by —1. If our program did not require the variables to be non-negative, we can 
always replace a variable a; by a;+ — a;_ , where x+ and x_ are both required to be 
non-negative. If the program had some inequalities of the form ai^'^x > a, or of the 
form ai^x < a we can first convert all the inequalities of the first type to those of 
the second type by negation, and then convert them to equations by introducing 
slack variables Xi > 0, and requiring aj^'^x-l-Xi = a. Similarly, any pr ogra m can be 



made to look like the dual program P* in the statement of Theorem 3.2, 

Program L is almost in the primal form needed by the duality theorem, but for 
two differences: there is no objective function, and we want the primal variables 
(the angles of the triangles) to be strictly positive, rather than just non-negative. 

For the moment, let us sweep these issues under the rug, by setting the objective 
function to be 0, and allowing the angles to vanish - it will be quickly apparent 
how to fix things up later. Let the modified program be Li. The dual LI of Li is 
the following: 

The dual program. Maximize F(u, v) : 7rJ2teF{T) + ^eeE{T) ^i^)'"e, subject 
to the conditions ut + < whenever e is an edge of t. 



Theorem 3.3. Assume that the conditions described in conditions and 2. 
are satisfied. Then the objective function of L\ is nonpositive. It equals if and 
only if there is a u, such that Ut = —Ve = u, for all t G F{T), e G E{T). 

Proof. Observation 1. Note that if there is a u as required in the statement 
of the Theorem, the objective function is, indeed, equal to 0. This is nothing other 



than the equality case of Condition 2.2 



Now, let u — min(Mi, . . . , Let u^P — Ui — u, and let v'^p = v 



u. 



for all values of the indices. The new variables are still feasible for L|, and by the 
observation 1 above, this transformation does not change the value of the objective. 
Furthermore, if ul^-* = for all i, then the objective is non-positive, and is equal to 
zero only if all of u^^'' are equal to zero as well - this is so, since all of the Vj must be 



non-positive, and all of their coefficients are positive, by Observation 2J. Assume, 
then, that u^^^ > for i G J^^^^; T^^^ is a proper subset of F{T) by construction. 

Observation 2. Suppose now that ~ u, for all t £ J-\. Then 

< TT ^ u- ^ A(e)u < 0, 
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by I 

Now, let u^^'^ — min(gjr(i) (wj"'^''), and let — u'l^ = if t ^ JF^^), and otherwise 
uf^ = ui^ — M^^). Likewise, vf^ = v"}'^ if e ^ E(T'^^^ , and otherwise v^^ = 

This still leaves us in the feasible region of L\ and strictly increases the value of 
the objective (by Observation 2). The new nonzero set JF^^^ is a proper subset of 
and we can repeat this process. In the end, we will wind up with a feasible 
point u^'^^jV*^'^^ with u^'^-' = 0, where the value of the objective is non-positive 
(by Observation 1), and strictly greater than the value of the objective at u, v (by 
Observation 2), thus completing the proof. □ 

The above theorem shows that the feasible region of L\ is non-empty. In order 
to find a solution with strictly positive angles, we write our angles ai as ai = (3i + e. 
We require all of the (3i to be non-negative, and our new objective is simply — e. 
Call the resulting program L2. Its dual L2 has the following form: 

Second dual. Maximize F(u, v) : ti" X]teF(T) + See£;(T) subject to 

the conditions ut + Ve < whenever e is an edge of t, and also 3^^^^^^^ ut + 

Theorem 3.4. The optimal value of the objective of L2 is strictly negative. 





Proof. Suppose the contrary. Then ut = —Ve = u, for some m, by Theorem |3. 3. 
However, in that case the last inequality of L2 is not satisfied, since the left hand 
side of the last constraint of vanishes. Indeed, it is equal to 

e(£E{T),e^dT 

Observe, however, that 3 X]tG-F(t) ^ counts each non-boundary edge with multiplic- 
ity 2, and each boundary edge with multiplicity 1. □ 

4. DIHEDRAL ANGLES OF DELAUNAY TRIANGULATIONS 

In the preceding section it was shown that in order for the linear program L to 
have a solution, it is necessary and sufficient that the putative dihedral angles 5{e) 



satisfy the inequaliti es and Below, we use these to derive another set of 
necessary conditions |4.3| , and show that under the additional assumption that the 



original dihedral angles do not exceed tt, these are equivalent to the inequalities 2.2 



The virtue of the inequalities |4.3| is that it is easier to interpret them geometrically. 



A special case of them is one of the main results of |20 - the connection will be 
explained in section |^. 

First, some definitions (these are the analogs of the definitions of section ^ without 
reference to face angles) - as before T is a triangulation of S, V{T) is the set of 
vertices of T, E{T) is the set of edges, and F{T) is the set of faces. We assume 
that each edge e G E{T) is given a weight (5(e) G R. The weight 5{e) will be called 
the dihedral angle at e. The quantity 6'{e) = n — (5(e) will be called the exterior 
dihedral angle at e. 
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If u is a vertex of T, then the cone angle Cy is defined to be 

e incident to v 
while the curvature Kiy IS defined to be Ky = 2t: — Cy 

Note. The cone angle of a boundary point is, thus, not the same as the boundary 
angle (in the language of section |l|, but smaller by tt. 

Below, it will often be useful to talk of the Poincare dual of T. Recall that this is 
the complex T* , such that the set of vertices of T* is in one-to-one correspondence 
with the set of faces of T, the set of edges of T* are in one to one correspondence 
with the edges of T - two vertices of T* are joined by an edges if and only if the 
corresponding faces of T share an edge, and finally the faces of T* correspond to the 
vertices of T - the vertices of a face v* of T* correspond to the faces of T incident 
to the corresponding vertex v. 

Definition 4.1. A subcomplex J- of T is closed, if whenever a cell t is in T, 
then so are all of the lower- dimensional cells incident to t. 

Definition 4.2. The total curvature of a subcomplex J- of T is defined as 

Notation. Let be a subcomplex of T, and let E'{!F) be the set of those edges 
e of T which are not edges of faces of T, but such that at least one endpoint of e 
belongs to T. For each edge e of T, define nj^{e) to be the number of endpoints of 
e which belong to J-. 

For example, if e e E{T), then njr(e) = 2; if e ^ E{T) U E'{T), then njr(e) = 0. 

Theorem 4.3. For every non-empty subcomplex J- ofT the following are equiv- 
alent: 

(a) . 

ee_E'(jr) 

with equality if and only if J- = T. 

(b) . Conditions [I^ and [E^ hold. 

Proof. We will show that (b) (a); the converse is immediate. 

vev{J^) vee vev{3^)v<Ee 
The last sum of equation |^ can be rewritten thus: 

Y E ^'(^)=2 E ^'(^)+ E Me)^'(e). (10) 

''ev'(.?-) e incident to w eGB(:F) eeE'{j=-) 

Finally, using the definition of S'{e), and combining equations ||and ^it follows 
that 

K{J^) = 2nx{T) + 2{ J2 S{e)~n\F{:F)\)- J2 '^'(e). (H) 
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Now, assume that the dihedral angles satisfy the inequalities 2.2. This means that 



the middle term on the right hand side of equation 11 is non- negative (and strictly 
positive unless T = T). □ 

4.1 Some corollaries and refinements 

Here are some easy consequences of Theorem ^7^ : 

Special case 1. !F = T. Then, since E'{!F) = 0, we just get the Gauss-Bonnet 
theorem (with curvature defined in terms of the dihedral angles). 

Special case 2. The complement of J- is an annulus containing no vertices of T 
- this corresponds to a simple cycle in the Poincare dual T* . For every edge in 
E'{T), Tie = 2. Theorem U just says that Y^eeE' ^'(e) > 0. 

Special case 3. For every edge in E'{J-), Ue = I. This will hold in the case where 
!F is a subset of T with a collar, and removing the edges in E' separates the 1- 



skeleton of T into (at least) two connected components. In this case. Theorem O 
tells us that 

J2 5'{e)>2nx{:F)-K{T). 



A drawback of both conditions 2.2 and Theorem 4.3 is that they require their 
respective inequalities to be checked for every subcomplex of T in order to verify 
whether a given assignment of dihedral angles is admissible. It is not hard to see 
that this requirement can be weakened somewhat. 

Let J- be the subcomplex in question. 

Observation 4.4. It can assumed that every edge in E{J-) is an edge of at least 
one face. 

Observation 4.5. The 1-skeleton of the Poincare dual J-* of J- can be assumed 
connected - this is somewhat stronger than saying that T is connected. 

Observation 4.6. It can he assumed that the 1-skeleton of J- - the Poincare 
dual of the complement of J- - has no isolated vertices. 

It is straightforward to check all of the above observations. 

In the special situation where (5(e) < tt for all e S E{T) - that is, S is Delaunay, 
one can further assume that every face / of .7^ is adjacent to at most one face of 
J^. Otherwise, adjoin / to to create a new complex !F' . This complex has one 
more face than !F, but its sum of dihedral angles is at most tt greater than that of 



J^. Hence, it is enough to check that satisfies the hypotheses of Theorem 4.3, or 



the conditions 2.2 



Definition 4.7. A simple subcomplex of T is a subcomplex T such that both T 
and T\J- are connected. 

Theorem 4.8. To verify that Conditions 2.1 and 2.2 (or equivalently, conditions 
of Theorem \(.^ (b)]) hold for all subcomplexes of T , it is necessary and sufficient 
to check them for simple subcomplexes. 
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Proof. By Observation 4.5, it is enough to check the connected subcomplexes of 
T. If such a subcomplex T is simple, then we are done. Otherwise, its complement 
is not connected. Let C be a connected component of the complement of T . 

Lemma 4.9. The complex C is simple. 

Proof. By construction, C is connected. Also, every point of T\C can be con- 
nected by a path to a point of J-. Since is assumed connected, the lemma follows. 
□ 

Consider T' — T\JC. By assumption, T' ^ T. Now 



^ ^{e)-i:\F{T)\ 



A(e)-7r|F(C)|l -J^C (12) 



A(e)-^|F(.F)'| 

,eG-E(.F') 



By lemma 2.5 (more precisely, a version for simple complexes). 



(13) 



Y A(e)~^|i^(C)|| ~EC<0' 



thus. 



Y A(e)-vr|F(^)| 



J \eeE(r') 



A(e) -'K\F{TYrime\ 



(14) 



(15) 



Thus, in order to check that condition p.2| holds for JF, it is enough to check 
that it holds for JF'. The proof of Theorem |4.8| is finished by the obvious induction 
argument on the number of connected components of the complement of J-. 

In the case where T is a genuine simplicial complex (that is, two cells intersect 
in a lower-dimensional cell), simple subcomplexes corresponds to non- coterminous 
minimal cutsets of edges of T: 

Definition 4.10. A collection C of edges of T is a cutset if removing the edges 
in C disconnects the 1-skeleton of T . A cutset C is minimal if no proper subset of 
C is a cutset. A cutset C is coterminous if all of the edges in C are incident to the 
same vertex. 

In other words, a minimal cutset corresponds to a separating simple curve in 
the Poincare dual T* of T (the curve need not be closed if T has boundary). A 
coterminous cutset corresponds to a boundary of a face in the dual. The non- 
coterminous simple cutset corresponding to the subcomplex J-' is nothing other 



than the set of edges E'{!F) define just before the statement of theorem 4.3 

In the case where the surface 5 is a flat disk, theorems 3.3, 1.6, 4.3, and [4.8 
immediately imply: 
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Theorem 4.11. Let T be a triangulation of the disk, let A : E{T) (0, tt] be 
an assignment of dihedral angles to the edges of T , let Vg (T) be the set of boundary 
vertices ofT and let A : Vq{T) —>■ (0, tt] be the assignment of boundary angles. Then, 
there exists a collection of points pi, ... ,pv(T) in the plane whose Delaunay 
triangulation is combinatorially equivalent to T , has dihedral angles given by A, 
and whose convex hull is a polygon with angles given by A, if and only if: 

— // V is an interior vertex of T , then 

(tt - A(e)) = 27r. (16) 



E 



e S E{T) incident to v 
-If V is a boundary vertex of T, then 

J2 (tt - A(e)) + (tt - A(w)) = 27r. (17) 

e G E{T) incident to v 



-If £ is a non- coterminous minimal cutset corresponding to a simple closed curve 
in the Poincare dual T* of T , then 

^(tt - A(e)) > 2^. (18) 

-// £ is a non- coterminous minimal cutset not corresponding to a simple closed 
curve in T* , and £ separates the boundary vertices ofT into two groups vi, . . . ,Vk 
and Vk+i, ■ . . ,v\Vo{T)\, then 

k 

^(7^-A(e))+^(^-A(^;0)>2^. (19) 



Remark. Theorem |4.1l| is nothing but one of the main resuhs (Theorem 0.1) of 
[^o| , stated in a different language (without mentioning convex ideal polyhcdra) , so 
we have succeeded in deducing that theorem in a purely combinatorial way from 
the results of jl^]. In addition, Theorem |4.8| is seen to be a direct generalization of 
pO| [Theorem 0.1] to a characterization of dihedral angles of Delaunay triangulations 
of arbitrary, possibly singular, Euclidean surfaces. 

5. A NETWORK FLOW APPROACH 



There is an alternative way to prove Theorem 3.3 which uses the Max Flow-Min Cut 
theorem of network flow instead of the duality theorem of linear programming. It 
should be noted that the difference between the two arguments is largely superficial. 



since the proof of Theorem 3.3 can be seen to essentially prove the Max Flow-Min 
Cut theorem. There are two reasons to set up the question as a result on network 
flow. The first is that the proof (hopefully) becomes clearer and more intuitive. 
The second is that the special sorts of linear programs that arise in the theory of 
network flow have been heavily analyzed from the viewpoint of complexity, which 
will allow us to give a very satisfactory estimate (section ||) for the running time 
of an algorithm to determine whether a structure with prescribed dihedral angles 
actually exists. 
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A network is a directed (multi)graph N, with two distinguished vertices - s (the 
source), and S (the sink). Each edge of N has a certain capacity, which is a real 
number, which is an upper bound on the amount of the commodity which can flow 
from the tail to the head of the edge. A cutset C of iV is a collection of edges, the 
removal of which leaves s and S in two different connected components of N\C. The 
capacity of the cutset C is simply the sum of the capacities of the edges comprising 
C. The capacity of an empty collection of edges is, of course, 0. 

The Max Flow-Min Cut theorem of network flow (see, for example, |2^, Chapter 
7]) says the following: 

Theorem 5.1 (Max Flow - Mm Cut). The maximal amount of a commod- 
ity that can flow from the source s to the sink S of a network N is equal to the 
capacity of the smallest cutset in N . 

This theorem can be proved in a number of ways. The interested reader can 



adapt the proof of The orem |3.3| to show Theorem 

To use Theorem 5.1 for our purposes, we need to set up a network of a special 
sort, starting with a triangulation T. 

The vertices of this network are divided into four classes: {s}, F{T), E{T), and 
{S}. 

The source s is connected to all of the vertices corresponding to the faces of T 
with edges of capacity 1. We call them edges of level 1. Every vertex corresponding 
to a face t is connected to the three vertices, corresponding to the three edges of 
t by edges of capacity 1. These are edges of level 2. Finally, each e G E{T) is 
connected to the sink S by an edge of capacity (5(e). These are edges of level 3. 

The following statement is self-evident: 

Observation 5.2. There exists a solution of the linear program Li if and only 
if the maximal flow through the above- constructed network Ni is equal to \F{T)\. 



Proof of theorem |3Jj. Consider a cut C of iVi. This will have some edges at 
level 1, removing which which will disconnect a subset Fq of the faces of T from the 
source. It is then not necessary to remove any edges of level 2 emanating from Fq. 
Let F{T)\Fo = Fi. Let F2 C Fi be those faces / for which the cutset contains all 
three edges of level 2 emanating from /. Finally, F3 = Fi\F2. All of the edges of 
level 3 (indirectly) emanating from F3 must be in the cutset C. These are precisely 
the edges corresponding to the edges of the subcomplex of T whose faces are in F3. 

What is the capacity of C? Evidently, it is equal to 

eeE(F3) 

If we want the flow through A'^i to be F, we must have 

Fo + 3F2+ ^(e) ^ ^' 

Or, noting that Fi = F - Fq, 
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In the special case where F2 = 0, it foUows that 

E ^(^)>P- (20) 

eG-B(-F3) 



Since F3 could have been any subset of F, it follows that condition 2.2 is necessary 
(this is, in any event, self evident). 



On the other hand, if condition 2.2 holds for any subcomplex of T, substituting 
the inequality 2.2 into |[ we see that 

c(C) >Fo + 3F2+F3^F + 2F2 > F. 

□ 

The above result is not quite what is required: the non-strict inequality ^ implies 
the existence of a consistent assignment of angles to the faces of T, but some of 
these angles may be equal to 0. In order to have the angles strictly positive, we 
must modify the weights in the network A^i as follows: 

Modify the capacity of level 1 edges to be 1 — 3e; those of level 2 to be 1 — e, 
and those of level 3 to be 1 — 2e. Call the resulting network N2. The existence 
of a consistent assignment of angles to the faces of T where all the angles are no 
smaller than e is obviously equivalent to the maximal flow in being equal to 
F(r)(l-3e). 

Consider now a cut C in A^2, with notation as before. The capacity of C will be: 
(l-3e)^^o + 3(l-e)F2+ ^ S{e) - 2eE{F3). (21) 

eeE{F3 



Assume that F2 = 0. The MinCut condition together with expression |2lj gives: 
(1 - 3e)Fo + E '5('^) - 2e^(^3) > (1 - 3e)F, 

eeE{F3 

or 

E ^(e) > (1 - 3e)^^3 + 2eE{Fs). 

e£E{F3 

li^cF^cF, then ^(^^3) > |i^3, and so EeeEiF^ ^i^) > F3. 

Suppose now that 6{F') — F' > ip > 0, for all proper non-empty subsets F' of F. 
Then the expression |2^ is no smaller than 

(1 - 3e)Fo + 3(1 - e)F2 + F3 + iP ~ 2eE{F3). (22) 

The edges of F3 can be divided into two classes: interior edges (those incident 
to two triangles of i^3) - these number Ei{F3), and boundary edges of F3 - those 
incident to only one triangle. These number Eq{F3). Since F3 is a proper subset 
of i^, EoiFs) > 0. Clearly, 

2E{F3) = 3F3 + Ea{F3). (23) 

The lower bound |2^ can thus be rewritten as 

c{C) > (1 - 3e)Fa + 3(1 - e)F2 + (1 - 3e)F3 + - eEgiFs) (24) 
- (1 - 3e)F + V + (1 - 3e)F2 - eEaiFs). (25) 
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Since e < ^ (a triangle cannot have all angles greater than 7r/3), we get 

c{C) - (1 - 3e)F >i,~ eEoiFa). (26) 

In other words, if J2e£E(F>) ^i^) ~ F' > ijj for every % <Z F' d F, there is a 
solution of the linear program Li with all face angles of all triangles no smaller 
than ip/E{T). □ 

6. DELAUNAY TRIANGULATIONS OF INFINITE SETS OF POINTS 

In this section we will show that theorem for singular Euclidean structures can 
be extended without change to infinite locally finite complexes: 

Theorem 6.1. LetT be an infinite but locally finite complex, and let A : E{T) —>■ 
(0, tt]. Then there exists a singular euclidean structure on T, with cone points at 
vertices ofT, whose Delaunay triangulation is combinatorially equivalent to T, and 
whose dihedral angles are given by A if and only if each finite subcomplex J- <Z T 
with E{J-) ^ has positive excess. 

There are two ingredients in the argument. The first (Lemma |6.2| ) is an extension 
of section |[ the second (|6.3| ) is a geometric estimate which will enable us to extract 
the necessary subsequences. 

Lemma 6.2. Let T be a complex, and let A : E(T) — > (0,7r]. Then there exists 
a Euclidean structure with cone angles at vertices of T, and dihedral angles given 
by A, except at the boundary edges of T, where they are smaller than prescribed by 
A if the excess of any subcomplex J- of T , such that E{T) ^ ^ is positive. 



Remark. Lemma 6.2 can be viewed as a relative version of theorem K.l 



Proof. The argument parallels very closely that of section ||. The existence of 
the desired structure is, as before, equivalent to a negative objective of a linear 
program, and as before, we set up a slightly simpler linear program first. To wit, 
the program L' is: 

Fositivity. All angles are strictly positive. 
Euclidean Faces. For every face t, at + (it + It = tt. 
Boundary edge dihedral angles. For every boundary edge e ~ AB, in- 
cident to the triangle ABC , 7 + = A(e), where the sla ck variables 



(see the comments following the statement of Theorem 3.2) Xg 0/rG also 
non-negative. 

Interior edge conditions. For every interior edge e = AB, incident to 
triangles ABC and ABD, a + 5 = A(e). 

We relax the program L' to a program L'l by dropping the requirement that the 
angles be strictly positive, make the objective 0, as before, and we see that the dual 
to the new weakened linear program L'l is the following: 

The dual program. Maximize F(u, v) : T^^i^pij-) ut + J2e£E{T) ^{^)'"e, subject 
to the conditions 

the inequalities of L^. + We < whenever e is an edge of t. 
new inequalities. Whenever e is a boundary edge of T, Ve < 
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Since the constraints of the above program V* are a superset of the constraints 



of L'l, Theorem 3.3 still tells us that the objective function is maximized if there 
exists a u, such that ut = —Ve = u, for all t G F{T), e € E{T). Now, this is not 
enough to guarantee that the objective is zero, since the equality case of condition 



2.2 (when J- — T) no longer exists. Since the new inequalities require u to be 
non-negative, it follows that for the objective function to equal 0, u must be 0. 

Now, we follow Section |3| again, to define a program L'2 in the same way as before 
(that is, since we want the angles to be strictly positive, we set ai — (3i + e, etc, 
and to also define its dual Lj*- By the same reasoning as before, it follows that the 
objective of L'2 , and hence of Lji is negative. In fact, we can do more: we can also 
require all of the slack variables Xe to be strictly positive. The (yet another) new 
dual program L3* will have the form: 

Third dual. Maximize -F(u, v) : tt X]te_F(T) + '^eeE(T) '-^(^)'^<^^ subject to the 
conditions 

Ut + Ve < whenever e is an edge of t. 
We < 0, when e is a boundary edge of T. 

3 Y,teF{t) "t + 2 J2eeE{T),e<fdT + 2 J2eeE{T) ,eedT < -1- 

If we omit the last constraint, we remain with the dual program L'j*, and as the 
discussion above showed, the objective of that can only be ii ut = Ve = 0, which 
is at odds with the last constraint of L3* □ 

Theorem 6.3. Let pi, . . . ,p„ be a set of points in the plane, and D their Delau- 
nay triangulation. Assume that the shortest edge of D has length 1, and the excess 
of every non-trivial subcomplex of D is no smaller than d^. Then 



diameter D < 



The proof of Theorem 6.3 will depend on a couple of easy auxiliary results: 



Lemma 6.4. Let ABC be a triangle with a/c,a/b < e, < e < 1/10. Then, 
a < 2e. 

Proof. This is an immediate consequence of the theorem of Cosines, or the 
theorem of Sines. □ 

Lemma 6.5. Let J- be a subcomplex of D. Then the excess of J-" is 

liABC), 



E 



triangles ABC such that AB G J- 
where j{ABC) is the angle at C. 

Proof. This is immediate from the definition of excess. □ 



Proof of theorem 6.3. Suppose that the conclusion of the theorem does not 
hold. Assume, without loss of generality, that the edge between the vertices po and 
pi is the shortest one (and thus of length 1). 



We construct a family of disks Di, D2, ■ ■ ■ , Dn, all centered on pq, and such that 

In \ 
d 0) 



the radius of Di is equal to (^n)"- ^» ~ L)i+i\Di, and let Ti to be the 
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maximal closed subcomplex of D contained in Di. The hypothesis of the theorem 
ensures that at least one of the annuli Ai contains no vertices of let this annulus 
be Aj. Then we claim that the excess of is smaller than do. Indeed, consider a 
triangle ABC of D adjacent to J-j along an edge AB. The vertex C of ABC lies 

outside Dj+i, and thus the lengths of AC and BC are at least (^^^ 
while the length of AB is at most (^^^ ■ Thus, 

2 _ do 

I n' 



do 



j{ABC) < < 



do 



by lemma 6.4. Thus, by lemma 6.5, it follows that 

excess J-j < do , 

contradicting the hypothesis of the theorem. □ 



Proof of Theorem 6.1. We only need to show that the positive excess condi- 
tions are sufficient, since they are obviously necessary. In addition, we may assume 
that the 1-skeleton of the Poincare dual T* is connected (if not, we prove the 
theorem for each connected component separately). 

Now, we pick a pair of adjacent base vertices vi,V2 £ V{T), and fix ^2) = 1. 
Now, for v,w & V{T) we define dc{v,w) to be equal to the combinatorial distance 
between v and w in the 1-skeleton of T. Now, define J^i to be the span of all 
vertices u, such that dc{vi, u) < i. The complex is finite by local finiteness of T. 
In addition, [J^^J^i = T, so every finite subcomplex of T belongs to some J^i. For 
each !Fi w e co nsider a geometric realization S{J-'i), whose existence is guaranteed 
by lemma |6.2| (there are many such realizations, we pick any one of them). 

Now, enumerate the faces of T, in such a way that ti contains the edge V1V2, and, 
for any j, the faces tj and tj+i are adjacent (that is, share an edge) in T. For each 
triangle, we have the space of shapes (similarity classes), given (for example) by the 
complex parameter z, obtained by placing the first two vertices of the triangle at the 
points and 1 in the complex plane, and reading off the position of the third point 



(in the upper half-plane, assuming the triangle is positively oriented). Theorem 6.3 
tells us that for any face t of T, the set of shape parameters of realizations of t is 
contained in a compact set Ct (since the ratio of lengths of any two sides is bounded 
by some constant, depending on the function A). We can think of each S{J^i) as 
being an element of C = Yij ^t, , which is a compact set by Tykhonov's theorem, 
and hence we can extract a convergent subsequence from S{Ti), . . . , Sf^Tk), ■ ■ ■■ Call 
the limit of that subsequence S. Since the dihedral angles are obviously continuous 
functions of the triangle parameters, the dihedral angles of S will be given by A, 
and so S is the sought-after realization. □ 



Note. For the comfort of more analytically inclined readers, it should be pointed out 
that the last argument is just an Arzela-Ascoli - uniform convergence on compact 
sets argument. 
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6.1 Remarks and Corollaries 

Let us assume that the complex T has no boundary, and that the function A is such 



that all of the cone angles (computed using Lemma 1.2) are equal to 27r, so that 
any realization is Euclidean (and hence can be developed into the plane). Then, 
we can combine theorems 3.1 and 4.11 to get 



Corollary 6.6. Let T he an infinite locally finite triangulation, and let A : 
E{T) —>■ (0, tt] be an assignment of dihedral angles to the edges ofT, let Then there 
exists aflat surface S, and a collection of points pi, . . . ,p„, . . . in S whose Delaunay 
triangulation is combinatorially equivalent to T and which has dihedral angles given 
by A if and only if: 



E 



-// V is an interior vertex of T , then 

(tt - A(e)) = 27r. (27) 
e G E{T) incident to v 

— If £ is a non- coterminous minimal cutset, then 

^(7r-A(e))>2^. (28) 

The above Corollary can be viewed as an extension of Theorem 0.1 of to 
the case of ideal polyhedra with infinitely vertices (and hence also of Andre'ev's 
theorem for ideal polyhedra Q), but not without certain caveats: even when T is 
topologically a disk, it is not at all obvious whether the metric on the surface S (or, 
even, any surface S) is geodesically complete. If it is complete, it follows that S is 
the Euclidean plane, and that the developing map is a global isometry, but otherwise 
the developing map is, only an immersion. Checking if the given simply-connected 
Euclidean surface S is the Euclidean plane is nothing other than the type problem 
- that is, we want to know whether a Riemann surface is parabolic, hyperbolic, or 
elliptic. If we know the shapes of all the triangles, this can be shown to be equivalent 
to the recurrence of a random walk on the 1-skeleton of the complex T, where an 
edge AB, incident to triangles ABC and ABD has weight l/(cotC-|-coti:)) (Q). 
Using this, it can be shown without too much difficulty that in the case where all 
of the dihedral angles are rational multiples of tt, bounded away from and tt, then 
this is equivalent to the recurrence of the symmetric random walk on the 1-skeleton 
of T. 



In the special case of all dihedral angles being equal to n or tt/2, corollary 6.6 
reduces to an existence theorem for infinite locally finite circle packings. In this 
case, the type problem, and a number of others, has been studied at great length 
by a number of authors, starting with Koebe, but more recently by A. Marden (in 
the context of Schottky groups), W. Thurston, B. Rodin and D. Sullivan, Z.-X. He, 
O. Schramm, and others. 



It should be noted the Theorem S.l says nothing about uniqueness, and the proof 



certainly does not show any form of uniqueness. In view of Theorem 1.4 one might 
suspect that perhaps this could be shown with more work. In fact, it is quite clear 
from the above-mentioned work on circle packing that uniqueness fails, though in 
a controlled manner described in the conjecture below: 
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Conjecture 6.7. LetT be a infinite locally finite complex, A a system of dihe- 



dral angles satisfying the hypotheses of Theorem 6.1. If there is a realization of T 
supported on a simply- connected domain Qq C C, then there is one supported on ev- 
ery simply- connected domain C C. Furthermore, such a realization is determined 
uniquely by fl (up to the group of Mobius transformations fixing fl). 

7. DELAUNAY CELLS IN MODULI SPACE 

In the Introduction we alluded to the cell decomposition of the moduli space A4 
of Euclidean structures, where the cells are given by Euclidean structures where 
the Delaunay triangulation has a fixed combinatorial type T. Each cell is a convex 
polytope P{T), and the results above can be used to describe the combinatorial 
and geometric structure of P(T) in some detail, although some interesting questions 
(discussed at the end of this section) remain open. 

Consider a codimension 1 face / of P{T). The face / may be either a boundary 
face of the moduli space M (viewed as a polyhedral complex) or an interior face. 
In the latter case, / corresponds to a change of combinatorial type of Delaunay 
triangulation from T to T', and it is well understood that the primitive such change 
is given by a diagonal flip, so / corresponds to a cyclic quadrilateral ABCD, where 
in T the Delaunay triangulation has triangles (for example) ABC and CD A, while 
in T', the triangles are DAB and BCD. On /, the quadrilateral ABCD can be 
triangulated either way, but the dihedral angle along (either) diagonal is equal to 

TT. 

Suppose now that / is a boundary face of - in particular, this will mean that 
none of the dihedral angles of the Euclidean structures in the interior of / are equal 



to TT. By Theorems 4.8 and 4.3, it follows that there is a simple cocycle c* of T* 
where the inequality 4.3[(b)] becomes an equality. This means, by the geometric 
estimates in the beginning of section ^, that the collar C of c* is becoming long and 
thin (that is, the conformal modulus of C diverges to oo), and so / corresponds to 
the Euclidean structure pinching off along the curve c* . The following construction 
(discussed in greater detail in an upcoming paper of the author) helps visualize this 
pinching off in terms of the more usual degeneration of hyperbolic surfaces: 

Consider a triangulated singular Euclidean surface (5, E, P). The data given by 
{S, E, P) (initially using a triangulation, though it is easy to show that the result is 
independent of triangulation) can be used to construct a cusped hyperbolic surface, 
as follows: 

(1) To each triangle t of the triangulation T of {S, E, P), associate an ideal triangle 
h{t). 

(2) For each pair of adjacent triangles ti — ABC and t2 — ABD of T, we have the 
log of the modulus of the cross-ratio of the four corresponding points: 

, , \AC\\BD\ 

(3) For each pair of adjacent triangles ti and t2 as above, glue the hyperbolic 
triangles h{ti) and h{t2) along the edge corresponding to AB with shear (see, 
eg. 111]) equal to r{ti,t2). 
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It is not hard to see that if we start with a Euchdean structure {S, E, P), we will 
wind up with a complete cusped hyperbolic structure (actually, it is sufficient, but 
not necessary, to start with a Euclidean structure - some similarity structures will 
give a complete structure also). The construction thus defines a map (certainly not 
injective, but which can be shown to be surjective using the construction of 0; 
|l^) between the moduli space of Euclidean structures with cone points and that of 
complete finite-area hyperbolic structures. It can be shown (this was an important 



part of 20 1 for the case of genus 0) that the degeneration, as above, of the 
Euclidean structure on (5, E, P), corresponds precisely to the pinching off along a 
simple closed curve of the hyperbolic structure on {S, H{E), P). 

Remark. Another surjection between the space of singular Euclidean structures 
and the space of cusped hyperbolic structures is well-known, and is discussed in the 
well-known paper of Troyanov p^ ]. The two surjections are not the same - this 
was remarked by C. T. McMuUen. 

8. COMPUTATIONAL COMPLEXITY 

Consider the following two decision problems: 

Problem 1. Let T be a simplicial complex, homeomorphic to a surface S (possibly 
with boundary), and let a : V{T) — > R"*" be an assignment of cone angles to the ver- 
tices of T. Does there exist a Euclidean structure E on S with the prescribed cone 
angles, such that the Delaunay triangulation of {S, E) is combinatorially equivalent 

to r? 

Problem 2. Let T be a simplicial complex, homeomorphic to a surface S (possibly 
with boundary), and let A : E{T) (0, tt] be an assignment of dihedral angles to 
the edges of T. Does there exist a singular Euclidean structure E on S, such that 
the Delaunay triangulation of {S, E) is combinatorially equivalent to T, and whose 
dihedral angles are given by A. 



By Theorem 1.6, there are efficient algorithms for both problems, since they 
reduce to the linear program L of section |l]. If the angles (cone angles in problem 1, 
dihedral angles in problem 2) are rational multiples of tt, such that the numerator 
and denominator are both bounded by C, then (by now) standard interior point 
methods allow us to solve the linear program L using Oiji^il + logC)) arithmetic 
operation, each involving arithmetic using precision 0(n(l -I- logC)). In the case 
where all of the prescribed cone angles are equal to 2-k, the logC can be disposed 
with, and we wind up with an algorithm of bit-complexity 0(n^ log^ n). 

Remark. In practice, the simplex algorithm appears much more efficient, and this 
has been used by M. B. Dillencourt to analyze all planar triangulations of up to 14 
vertices, and to determine which of them are combinatorially equivalent to planar 
Delaunay tessellations (|^; ||]). 

For problem 2, the network flow formulation of section p| turns out to give a 
markedly superior complexity. Indeed, it has been shown inllj that for a network 
with n nodes, m arcs, and the (integer) capacity of each arc bounded by U , we can 
determine the maximal flow in time bounded by 0(nm log((n/m)(log J7)^/^ -I- 2)). 

For the network A^i of section assuming the genus of the surface is fixed, 
the number of arcs and nodes in the network are both bounded by constant mul- 
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tiples of the number F of faces in the complex T. If aU of the dihedral an- 
gles are rational multiples of tt, with numerators and denominators all bounded 
in absolute value by C, the quantity U is bounded by C"^'^-' (since we need to 
compute the least common multiple of the denominators), giving a running time 
bound of 0(F5/2). p^j. ^^le program N2, the bound is the same, since the only 
difficulty consists of picking the right value of e, and this can be made to be 
1/lcast common denominator of the dihedral angles. 

9. SOME APPLICATIONS TO COMBINATORIAL GEOMETRY 

A well-known theorem of Steinitz says that every three-connected planar graph 
can be realized as the 1-skeleton of a convex polyhedron in i?'^, while a famous 
theorem of Aleksandrov states that every Euclidean metric on with positively 
curved cone-points can be realized as the induced metric on the surface of a convex 
polyhedron in R^. Below, we show a negative result, which should be compared 
with the final example of ||l7|| : 

Theorem 9.1. There exist infinitely many triangulations of S'^ which can not be 
realized as a Delaunay triangulation with respect to the cone-points of any Euclidean 
metric on S"^ with positively curved cone-points. 

First a definition: 

Definition 9.2. Let T be a triangulation. The stellation s(T) ofT is the com- 
plex obtained by replacing each face ABC of T by three faces AOB, AOC, and 
HOC. 



Theorem 9.1 follows immediately from the following claim. 

Claim 9.3. Let T be any triangulation of with at least eight faces. Then the 
stellation s{T) ofT is not combinatorially equivalent to the Delaunay triangulation 
of any Euclidean metric on S'^ with positively curved cone-points, where the cone- 
points correspond to vertices of s{T). 

Proof. Let an old vertex of s{T) be one that was already a vertex of T, while a 
new vertex be one that was added at stellation. The set M of new vertices of s(T) 
corresponds to the set of faces of T. For any vertex v, recall that C{v) denotes the 
cone angle at v. The Gauss-Bonnet theorem tells us that 

{2tt - C{v)) = 4tt. (29) 

veV{s{T)) 

Or, recombining the terms: 

Civ)^2ni\V{siT))\-2). (30) 

veV{s{T)) 

Note that every edge of s{T) is incident to an old vertex, and to at most one new 



vertex. Combining this observation with Lemma 1.2, we see that for any Delaunay 
triangulation combinatorially equivalent to s(T), it must be true that 

E ^(^)^ E ^(")- (31) 

vG{o\d vertices of s{T)} ti6{new vertices of s{T)} 
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Combining Eq. |T| with Eq. it follows that 

E Civ)>n{\VisiT))~2). (32) 

tie {old vertices of s{T)} 

Note now that \V{s{T))\ = \V{T)\ + \F{T)\. By a standard computation using 
Euler's formula for triangulations of the sphere, we know that |V^(r)| = ^\F{T)\ + 2, 
thus \F{T) \ = 2\V{T) \ - 4. Thus, 

\V{s(T)) -2\=3\V{T)\-6. (33) 

Equations |3^ and ^ together imply that the average cone angle at an old vertex 
of s{T) is at least 7r(3|y(r)| - 6)/\V{T)\ = in - 6/\V{T)\. If \V{T)\ > 6 it follows 
that the average cone angle at an old vertex is greater than 27r, which contradicts 
the assumption that the cone angles were positively curved. □ 

10. LINEAR HYPERBOLIC STRUCTURES ON 3-MANIFOLDS 

As explained in , the study of Euclidean triangulations on surfaces is essentially 
equivalent to the study of hyperbolic ideally triangulated complexes, which are 
combinatorially just cones over the triangulated surface. Hence, the contents of 
this section are closely related to the subject-matter of much of the rest of the 
paper, in more than just the linear programming approach. 

Consider a 3-manifold with boundary a collection of tori, and consider a 
topological ideal triangulation T of M^. We would like to know when there is a 
complete hyperbolic structure on M^, such that T is a geometric ideal triangula- 
tion. In general, this is a very difficult question, at least as hard as Thurston's 
hyperbolization conjecture (since even if admits a hyperbolic structure of finite 
volume, there might not be an ideal triangulation combinatorially equivalent to T). 
However, below we will consider a "linear" version of the question above. 

Recall that an ideal simplex S in has the properties that 

Euclidean links. The link of each vertex is a Euclidean triangle. 

Equal opposite dihedral angle. If = ABCD, then the dihedral angles at the 
edges AB and CD are equal (this is actually a consequence of the condition on the 
links). 

An ideal simplex is thus determined by the angles of the link of any one of its 
vertices (all links are easily seen to be the same). 

Now, if T comes from a genuine hyperbolic structure, it must be true that the 
sum of the dihedral angles incident to the edges of T must equal 27r, and so for T 
to correspond to such a structure, the following linear program must have a strictly 
positive solution: 

The variables. These are the dihedral angles of the simplices. For each simplex 
S we use three angles a, (3, 7 corresponding to the angles of the link of one vertices 
of S*. 

Simplex conditions. For each simplex 5, the sum of the dihedral angles is equal 
to tt: 



/s : a + /9 + 7 = TT. 



(34) 
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Edge conditions. For each edge e of T the sum of the dihedral angles of all sim- 
plices incident to e equals 2tt : 



Definition 10.1. We say that if the above linear program has a non-negative 
solution, then the pair {M^ ,T) admits a weak linear hyperbolic structure. // the 
linear program has a strictly positive solution, then we say that the pair (M^,T) 
admits a linear hyperbolic structure. 

It is clear that the existence of a linear hyperbolic structure is, in general, no guar- 
antee of the existence of a genuine hyperbolic structure, since the linear conditions 
do not preclude "Dehn surgery" singularities, as well as translational singularities 
along edges of T (see [|5|; |l^ for discussion). Conversely, as remarked above, the 
existence of a complete hyperbolic structure on is no guarantee that there is a 
positively oriented ideal triangulation combinatorially equivalent to T (or, indeed, 
any positively oriented triangulation). However, linear hyperbolic structures have 
the advantage of being considerably more tractble. 

Theorem 10.2. In order for there to exist a weak linear hyperbolic structure 
for {AI^,T), every normal surface with respect to T must have non-negative Euler 
characteristic. In order for there to exist a linear hyperbolic structure for {M^,T), 
every non-boundary-parallel normal surface with respect to T must have strictly 
negative Euler characteristic. 

Proof. We will use the method of Section ^. First, let us write the dual program 
of the linear program (referred to as Lh in the sequel) for weak hyperbolic structure: 
Our variables are {ws}, where S ranges over all the simplices of T, and {we} where 
S ranges over all the edges of T. 

The dual program LJj is; 

■ Maximize EseS(T) + 2J2e<£E{T) "e- 

. Subject to vs + + < for all faces S and pairs of opposite edges ei, 62- 

In order for the primal program to have a non-empty feasible region, the objective 
of the dual must be non-positive. 

Consider now a normal surface S (see, eg, [|[ |l^; ^ for rudiments of normal 
surface theory). The surface S intersects each simplex 5 in a collection of disks, 
which are combinatorially either triangles (cutting off one vertex of S from the 
other three), or quadrilaterals (separating one pair of vertices from another). For 
each simplex S, define ts{S) to be the number of triangular components oi S H S, 
and qs{S) to be the number of quadrilateral components. For each edge e of T, 
define ie{S) to be the number of intersections of S with e. The intersections of S 
with the simplices of T induces a triangulation T of S, where each triangular disk 
contributes one triangle, and each quadrangle contributes two. Define us{S) to be 
the number of triangles of T sitting inside a simplex S. Evidently, us = ts -\- 2qs. 

Note that, by Euler's formula, x('5) = J2eeE{T] — \ SseS(T) '^S- This is seen 
to be very similar in form to the objective function of Zj^, so let us set vs — —us, 
and Vf, = ie- 




(35) 
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Lemma 10.3. The assignment of the variables as above satisfies the inequality 
constraints of L*^ . 

Proof of lemma. We need to check that for S a simplex and 61,62 a pair of 
disjoint edges of S. We need to check that 

iei + ie2 < Us- (36) 

By linearity, we need just check the inequaUty ^ for connected components oiSOS. 
If that component is a triangle t, then t contributes 1 to ie^ + (since a "normal 
triangle" intersects exactly one of each pair of opposite edge). Also, t contributes 
1 to Us, so for a triangular face, the right and left hand sides of (|3^) are equal. 

Suppose now we have a quadrilateral component q. It contributes 2 to the right 
hand side of (|3^) . As for the left hand sides, q hits two pairs of opposite sides of 
S, so if 61 and 62 is one of those pairs, then we have a contribution of 2 to the left 
hand side, and otherwise we have a contribution of 0. □ 

Remark 10.4. Notice that if S is such that all of the components of S D S, 
for all S e S{T) are triangles, then all of the constraints of L'^ are equalities with 
the assignment of variables as above. Any such S is easily seen to be a union of 
boundary tori. 



Lemma 10.3 concludes the proof of the "weak" part of Theorem 10.2 , since if any 
S had positive Euler characteristic, the program L'^ would have a positive objective, 
and thus the program would have no solution. 

For the proof of the "strong part" we use the same trick as in Section ||. Define 
new variables a' — a + e, etc. Our primal linear hyperbolicity program Lg is now: 

. Minimize: — e 

. subject to face constraints a' + f3' + •y' + 3e ~ tt. 

. and to edge constraints J2aPrime +v{e)e — 2tt, where v{e) is the valence of e. 

The dual program L* is then 

■ Maximize Ese5(T) «s + 2Eee-E(T) "e- 

. Subject to the old constraints vs + + v,,^ < for all faces S and pairs of 
opposite edges 61,62. 

. and the new constraint 3j2sesiT) "-s + J2eeE(T) ^(6)^-E < — 1- 

In order for {M'^,T) to be linearly hyperbolic, the objective must be strictly 
negative. 

Observe that the sum of the left hand sides of the old constraints is equal to the 
left hand side of the new constraint. Indeed, each vs occurs three times (once for 
each pair of opposite edges), and each Ve occurs the number of times equal to the 
valence of e. Hence, the new constraints simply says that in at least one of the 



old constraints the inequality must be strict. Keeping in mind Remark 10.4, this 



implies that every non-boundary-parallel normal surface must have strictly negative 



Euler characteristic, thus proving the second part of Theorem 10.2 



Some remarks may be in order: It is easy to see (and not surprising) that an 
identical theorem can be proved if the cone angles around the edges of T are required 
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to not be smaller than 27r, while if the angles are smaller than 27r, one can show 
an analogous "orbifold" version of the theorem. A more interesting question is 



whether theconverse of Theorem L0.2 holds. This is equivalent to asking whether 
every assignment of variables satisfying the constraints of programs L'^ and L* 
comes from a normal surface. 
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